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Abstract This paper focuses on the branching process for solving any 
constraint satisfaction problem (CSP). A parametrised schema is pro- 
posed that (with suitable instantiations of the parameters) can solve 
CSP's on both finite and infinite domains. The paper presents a formal 
specification of the schema and a statement of a number of interesting 
properties that, subject to certain conditions, are satisfied by any in- 
stances of the schema. It is also shown that the operational procedures 
of many constraint systems (including cooperative systems) satisfy these 
conditions. 

Moreover, the schema is also used to solve the same CSP in different 
ways by means of different instantiations of its parameters. 
Keywords: constraint solving, filtering, branching. 



1 Introduction 

To solve a constraint satisfaction problem (CSP), we need to find an assignment 
of values to the variables such that all constraints are satisfied. A CSP can have 
many solutions; usually either any one or all of the solutions must be found. 
However, sometimes, because of the cost of finding all solutions, partial CSP's 
are used where the aim is just to find the best solution within fixed resource 
bounds. An example of a partial CSP is a constraint optimisation problem (COP) 
that assigns a value to each solution and tries to find an optimal solution (with 
respect to these values) within a given time frame. 

A common method for solving CSP's is to apply filtering algorithms (also 
called arc consistency algorithms or propagation algorithms) that remove incon- 
sistent values from the initial domain of the variables that cannot be part of any 
solution. The results are propagated through the whole constraint set and the 
process is repeated until a stable set is obtained. However, filtering algorithms 
are, often, incomplete in the sense that they are not adequate for solving a CSP 
and, as consequence, it is necessary to employ some additional strategy called 
constraint branching that divides the variable domains and then continues with 
the propagation on each branch independently. 
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Constraint Solving algorithms have received intense study from many re- 
searchers, although the focus has been on developing new and more efficient 
methods to solve classical CSP's [?,?] and partial CSP's [?,?]. See [?,?,?,?] for 
more information on constraint solving algorithms and [?,?] for selected com- 
parisons. 

To our knowledge, despite the fact that it is well known that branching step 
is a crucial process in complete constraint solving, papers concerned with the 
general principles of constraint solving algorithms have mainly focused on the 
filtering step [?,?,?]. 

In this paper, we propose a schema for constraint solving for both classical 
and partial CSP's that includes a generic formulation of the branching process. 
(This schema may be viewed as a generalisation and extension of the interval 
lattice-based constraint-solving framework in [?].) The schema can be used for 
most existing constraint domains (finite or continuous) and, as for the framework 
in [?], is also applicable to multiple domains and cooperative systems. We will 
show that the operational procedures of many constraint systems (including 
cooperative systems) satisfy these conditions. 

The paper is organised as follows. Section 2 shows the basic notions used 
in the paper and Section 3 describes the main functions involved in constraint 
solving with special attention to those involved in the branching step. In Sec- 
tion 4 a generic schema for classical constraint solving is developed and its main 
properties are declared. Then, Section 5 extends the original schema for partial 
constraint solving and more properties are declared. Section 6 shows several in- 
stances of the schema to solve both different CSP's and different solvings for the 
same CSP. Section 7 contains concluding remarks. Proofs of the properties are 
found in the Appendixes. 

2 Basic concepts 

Let -D, D\, . . . , D n be sets or domains. Then #D denotes the cardinality of D, 
p(D) its power set and D < denote any totally ordered domain. _l_£> and 
denote respectively, if they exist, the bottom and top clement of D and fictitious 
bottom and top elements otherwise. Throughout the paper, A denotes a set of 
domains called computation domains. 

Definition 1. ( Constraint satisfaction problem) A Constraint satisfaction prob- 
lem (CSP) is a tuple (V,T>,C) where 

— V = {vi, . . . , v n } is a non-empty finite set of variables. 

— V = p{Di) x ... x p(D n ) where D { G A. 

— C C p(Di, . . . , D n ) is a set of constraints for V. 

If, as in the above definition, T> = p(D\) x . . . x p(D n ), where Di £ A for 
alH e {1, . . . ,n}, then the set of all constraints for V is denoted as C-p and the 
set {Di | 1 < i < n} is denoted as Aj>. 



3 



Definition 2. (Constraint store) Let S = (di, . . . ,d n ) G V. Then S is called 
a constraint store for (V, T>,C). S is consistent if, for all i G {1, . . . , n}, di ^ 0. 
S is divisible if S is consistent and for some i G {1, . . . ,n}, jfdi > 1. Let 
S' = {d' x , . . . , d' n ) be another constraint store for (V, V, C). Then S ^ s S' if and 
only if di C d\ for 1 < i < n. 

S is a solution for (V, V,C) if S = ({si}, . . . , {s n }) and (s\, . . . , s n ) G c, for 
all c G C. S' is a partial solution for (V 7 T>,C) if there exists a solution S" for 
(V 7 T>,C) such that S" ~< s S' . In this case we say that S' covers S" . 

The set of all solutions for (V,V,C) is denoted as Sol((V,V,C}). Note that, 
if S G Sol((V, T>, C)), then 5* is consistent and not divisible. If (d\, . . . , d n ) G V 
and i G {1 . . . ,n}, then (di, . . . , d n )[di/d'\ = (d u . . . , d t -i,d', d i+u . . . , d n ). 

Example 1. Let V = p(Bool) x p(Bool) x p(Bool). Let c = x V (y A z) be a 
constraint for V. Then c = {(0, 1, 1), (1, 1, 1), (1, 0, 1), (1, 1, 0), (1, 0, 0)}. Let 

S 1 = ({1}, {0}, {0}), 5 2 = ({0, 1}, {0, 1}, {0, 1}), 

S 3 = ({0},{0},{0}), 5 4 = (0,0,{O}). 

Then, Si is a solution but 5*2, S3 and 5*4 are not. Note also that Si, S2 and S3 
are consistent and S4 is inconsistent. 

Definition 3. (Stacks) Let P = (Si,... , Sn) G p(T>). Then P is a stack for 
(V,V,C). 

Let P' = (S[,... ,S' e ,) be another stack for (V,V,C). Then P < p P' if and 
only if for all Si G P (1 < i < I), there exists S'j G P' (1 < j < £' ) such that 
Si dis S'j. In this case we say that P' covers P. 

3 The Branching Process 

This section describes the main functions used in the branching process. 

First we define a filtering function which removes inconsistent values from 
the domains of a constraint store. 

Definition 4. (Filtering function) filtering^ :: p(Cp) x T> — > V is a called a 
filtering function for V if, for all SeD, 

(a) filtering v (C, S) ;< s S; 

(b) VR G Sol({V,V,C}) : R^ s S =^ R ^ s filtering V (C, S). 

(c) If filtering V (C,S) is consistent and not divisible then filtering ^(C, S) is a 
solution for (V, T>, C) . 

Condition (a) ensures that the filtering never gains values, condition (b) 
guarantees that no solution covered by a constraint store is lost in the filtering 
process and condition (c) guarantees the correctness of the filtering function. 

Variable ordering is an important step in constraint branching. We define a 
selecting function which provides a schematic heuristic for variable ordering. 
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Definition 5. (Selecting function) Let S = (d\, . . . , d n ) E V. Then 

choose :: {S G V \ S is divisible] -> {p(-D) | D G Ap} 

is called a selecting function for V if choose(S) — dj where 1 < j < n and 
#dj > 1. 

Example 2. Here is a naive strategy to select the left-most divisible domain. 
Precondition : {S = (di, . . . , d n ) G V is divisible} 

choose nmve(S) = d 

Postcondition : {3j G {1, . . . , n} . d = dj , ffdj > 1 and 
Vi€{l,... ,j-l}:#di = l}. 

In the process of branching, some computation domain has to be partitioned, 
in two or more parts, in order to introduce a choice point. We define a splitting 
function which provides a heuristic for value ordering. 

Definition 6. (Splitting function) Let D e A and k > 1. Then 

split D :: p{D) -» p(D) x . . . x p{D) 

S v ' 

k times 

is called a splitting function for D if for all d e p(D), j^d > 1, this function is 
defined split D (d) = (d\, . . . ,dk) such that the following properties hold: 

Completeness : d\ U . . . U dk = d. 
Contractance : di C d, Vi G {1, . . . , k}. 

To guarantee termination, even on continuous domains, an extension of the 
concept of precision map shown in [?] is applied here. 

Definition 7. (Precision map) Let 5RX = (SR + , Integer) where 5R + is the do- 
main of non-negative reals. Then precision D is a precision map for D G A, if 
precision D is a strict monotonic function from p(D) to SRZ. 

Let S = (di , . . . , d n ) be a constraint store for (V, T>, C) and, for each D G A-p, 
precision D is defined for D. Then, a precision map for V = {D\, . . . ,D n ) is 
defined as 

precision(S) = ^ precision D .(di), 

l<i<n 

where the sum in is defined as (ai, 02) + (61, ^2) = (ai + «2 + ^2)- 

The monotonicity of the precision is a direct consequence of the definition 1 . 

1 5RX is continuous and infinite so that it is supposed that we can define a map from 
D to IKT, even if D is infinite. !KT was chosen since it is valid for the interval domain 
as it is shown in Section 6, but any domain totally ordered supporting the operator 
— may also be adequate (see Line 3 in Figure 1). 
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Proposition 1. Let S,S' be two constraint stores for (V,T>,C). If S -< s S' then 
precision(S) <ri precision(S'). 

The precision map also means a novel way to normalise the selecting functions 
when the constraint system supports multiple domains. 

Example 3. The well known first fail principle chooses the variable constrained 
with the smallest domain. For multiple domain constraint systems to emulate 
the first fail principle, we define choose /l so that it selects the domain with the 
smallest precision 2 . We denote this procedure by choose g. 

Precondition : {S = (d\, . . . , d n ) € T> is divisible} 

choose ff(S) = d 
Postcondition : {3j G {1, . . . , n) . d = dj , fj^dj > 1 and 

Vi € {1, . . . ,n}\{j} : #di > 1 precision D .(dj) <srj precision D .(di)}. 

4 Branching in Constraint Solving 

Figure 1 shows a generic schema for solving any CSP (V,T>, C). This schema 
requires the following parameters: C, the set of constraints to solve, a constraint 
store S for (V,£>,C), a bound p 6 and a non-negative real bound e. There 
are a number of values and subsidiary procedures that are assumed to be defined 
externally to the main branch procedure: 

— a filtering function filtering^/2 for T>\ 

— a selecting function choose /l for £>; 

— a splitting function split D for each domain D e Z\x>; 

— a precision map for V (therefore it is assumed that there is defined one 
precision map for each D 6 Zip); 

— a stack P e p(V) for (V,V,C). 

It is assumed that all the external procedures have an implementation that 
terminates for all possible values. 

Theorem 1. (Properties of the branch/4 schema) Let S be the top element in 
V (i.e., S = (Di, ...,D n )), e e 5R + andp = T^x- Then, the following properties 
are guaranteed: 

1. Termination: if £ > 0.0 then branch{C , S ', p , e) terminates; 

2. Completeness: if £ — 0.0 and the execution of branch (C, S, p, e) terminates, 
then the final state for the stack P contains all the solutions for (V,V,C); 

3. Approximate completeness: if e > 0.0 and R is a solution for (V, T>, C), then 
an execution of branch (C,S,p,e) will result in P containing either R or a 
partial solution R' that covers R. 

2 It is straightforward to include more conditions e.g., if di,dk, dj have the same (min- 
imum) precision, the most left domain can be chosen i.e., d minimum ( it k,j) ■ 
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procedure branch (C, S, p,e) 



S <- filtering ^{C, S); (1) 

if 5 is consistent then (2) 
if (S is not divisible or p < Tsrz and p — precision(S) < (e,0)) then (3) 

push(P, S); %% Add S to top of P (4) 

else (5) 

dj <— choose(S); (6) 

(dji, . . . , djfe) <— split D . (dj), where dj C IX,; (7) 

branch(C, S[dj/dji], precision(S), e) V ^ 

V \ %% Choice Points (8) 

branchiC, S[dj /djk], precision(S) , e) ; J 

endi/; 



end. 



Figurel. branch/4: A Generic Schema for Constraint Solving 



4- Correctness: if e — 0.0, the stack P is initially empty and the execution 
of branch(C,S,p,e) terminates with R in the final state of P, then R is a 
solution for (V, V, C) . 

5. Approximate correctness or control on the result precision: If Pq.o, Pa an d 
P £2 are stacks resulting from any terminating execution of branch(C, S, p, s) 
(where initially P is empty) when e has the values 0.0, e\ and S2, respectively, 
0.0 < £\ < £2 and Po.o is not empty, then Po.o dip Ps 1 dp Pe 2 ■ 
(In other words, the set of (possibly partial) solutions in the final state of the 
stack is dependent on the value of e in the sense that lower £ is, closer to 
the real set of solutions is). 

Observe that the bound e guarantees termination and allows to control the 
precision of the results. 

5 Solving optimisation problems 

The schema in Figure 1 can be adapted to solve COPs by means of three new 
subsidiary functions. 

Definition 8. (Subsidiary functions and values) let D < be a totally ordered 
domain 3 . Then we define 

— a cost function, fcost : : V — > D < ; 
3 Normally D< would be 3?. 
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— an ordering relation, o :: D < x D < G {>,<,=}; 

- a bound,, S G D < . 

Then the extended schema, branch+/A, is obtained from the schema branch /4 
by replacing Line 4 in Figure 1 with: 

if fcost(S) o S then (5 <- fcost(S); push(P, S) endif; (4*) 

Theorem 2. (Properties of the branch+fA schema) Let S be the top element in 
V (i.e., S — (Di, D n ) ), e G 5R + and p = Tsrx- Then, the following properties 
hold: 

1. Termination: if e > 0.0, then the execution of branch + (C, S,p,e) terminates; 

2. If f cost is a constant function with value 5 and o is =, then all properties 
shown in Theorem 1 hold for the execution of branch + (C, S,p,e). 

3. Soundness on optimisation; if e — 0.0, o is > (resp. <), 5 = -Ld < (resp. 
Td < ), the stack P is initially empty and the execution of branch + (C,S,p,e) 
terminates with P non-empty, then the top element of P is the first solution 
found that maximises (resp. minimises) the cost function. 

Unfortunately, if e is higher than 0.0, we cannot guarantee that the top of the 
stack contains a solution or even a partial solution for the optimisation problem. 
However, by imposing a monotonicity condition on the cost function fcost/1, we 
can compare solutions. 

Theorem 3. (More properties on optimisation) Suppose that, for i G {1,2}, 
P £i is a stack resulting from the execution of branch + (C, S, p, Si) where £j G 3? + . 
Suppose also that top(P) returns the top element of a non empty stack P. Then, 
if Si < 82 the following property hold. 

Approximate soundness: If for i G {1,2}, P £i is not empty, and top(P £2 ) is 
a solution or covers a solution for (V,T>,C), then, if fcost/1 is monotone and o 
is < (i.e., a minimisation problem), 

fcost(top(P Sl )) < D< f cost (top (P e2 )), 
and, if fcost/1 is anti-monotone and o is > (i.e., a maximisation problem), 

fcost(top(P £l )) h D< fcost(top(P £2 )). 

Therefore, by using a(n) (anti-)monotone cost function, the lower e is, the 
better the (probable) solution is. Moreover, decreasing e is a means to discard 
approximate solutions. For instance, in a minimisation problem, if 

fcost(top(P £l )) >- D< fcost(top(P £2 )) 

with fcost/1 monotone, then, by the approximate soundness property it is de- 
duced that top(P £2 ) cannot be a solution or cover a solution. 
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6 Examples 

To illustrate the schemas branch/^ and branch + /A presented in the previous two 
sections, several instances of branch /4 are given for some well-known domains 
of computation. In addition, we explain how the choice of instantiation of the 
additional global functions and parameters in the definition of branch + /4 can 
determine the method of solution for the CSP. 



6.1 Some instances 

In the following, branch x denotes an instance of the schema branch /4 for solving 
the CSP (V,V,C) where X C A v . We assume that 

A —{Bool, Integer, 5ft, Set Integer} U {Interv(D) D is a lattice}. 

where Interv(D) denotes the set {(di,d 2 ) | d\,d 2 € D,d\ < d 2 }. 

To identify branchy, we indicate a possible definition for both the splitting 
function and the precision map for each D e Ad and assume that both a se- 
lecting function (e.g., choose g as defined in Example 3) and a filtering function 
for V have been already defined. We also indicate the initial value of S G V, so 
that the execution of branchx>(C,S,p,s) allows to solve the CSP where e E 5ft + . 



The finite domain (FD) Constraint solving in a FD of sparse elements is 
solved by an instance branchFD as defined below where split FD is defined as a 
naive enumeration strategy in which values are chosen from left to right. For 
example, consider a finite domain of integers [?], Booleans [?] or finite sets of 
integers [?]). 



FD E {Integer, Bool, Set Integer}, 
S = (FD, . . . , FD); 



branchFD 



precision FD (d) = 0); 

split FD ({a 1 ,a 2 ,a 3 , . . . ,a k }) = ({a 1 },{a 2 ,a 3 , . . . ,a k }). 



Finite closed intervals Many existing FD constraint systems solve constraints 
defined in the domain of closed intervals [a, b] where a, b 6 FD and denoted 
here by a..b. Usually a,b are either integers [?], Booleans 4 [?] or finite sets of 
integers [?]. Here are two instances of our schema that solve CSP's on these 
domains: 

FD G {Integer, Bool}, 

S = (Lpd..~T pd, • • • , J-fd-.T fd); 



branch Inte rv(FD) < 



precision Interv{FD) (a..b) = (b- a,0); 

. SpUt Interv(Integer)( a -- b ) = (a. .a, a + 1 . .b) . 



The Boolean domain is considered as the integer subset {0, 1}. 
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FD = Set Int, 



branch Inte rv(FD) 



S = (0.. Integer, . . . , 0. .Integer); 

" v ' 

n times 

precision Interv(FD) (a..b) = (#6- #a,0); 
k split interv(FD)("-fy = (a~b\{c}, a U {c}..b) where c G 6\a. 



Lattice (interval) domain In [?], we have described a generic filtering algo- 
rithm that propagates interval constraints on any domain L with lattice struc- 
ture subject to the condition that a function o L ::LxL^5Ris defined that is 
strictly monotonic on its first argument and strictly anti-monotonic on its second 
argument. Below we provide an instance to solve any CSP defined on Interv(L): 



branch 



' L is a lattice and S = ([± L ,T L ], ... , [± L ,T L ]): 



Interv(L) 



precision Interv{L) {r) = { 

, SpUt Interv(L) 

({a,b}) 



n times 

{ (bo L a,2) ifr= [a,b]; 
(b o L a, 1) if r = (a, 6]; 
(bo L a, 1) if r = [a,b); 
k (b o L a,0) if r = (a, 6); 
({a, c], (c, 6}) where a c d. 



{a, 6} denotes any interval in L. With this instance we have a constraint solving 
mechanism for solving (interval) constraints defined on any domain with lattice 
structure. Thus it is a good complement to the filtering algorithm in [?]. Note 
also that if L is 5R and o L is — , we obtain the instance branch i nterv (jR) (also, if 
c = we have a usual strategy of real interval division at the mid point). 



A cooperative domain The schema also supports cooperative instances that 
solve CSP's defined on multiple domains. This is done by mixing together several 
instances of the schema branch /A. As an example, consider branch bnr as de- 
fined below where split ] nterv (D) an d precision i nterv (D) are defined as in previous 
examples for D G {Bool, Integer, 



(A = {Interv(D) | D G {Bool, Integer, ft}}. 
branch BNR j g = ( 0..p 1; ,„ ,0..£>» ), {D u ...,£>„} C {Bool, Integer, 5R}. 

This instance simulates the well known splitsolve method of CLP (BNR) [?]. 
The generic schema is also valid for solving non-linear constraints provided 
the filtering function filtering p/2 solves constraints in non- linear form. 
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6.2 Different ways to solve the instances of a CSP 

Here we show that, for any instance, the schema branch+/4 also allows a CSP 
to be solved in many different ways, depending on the values for fcost, S and o. 

For instance, a successful result for a classical CSP can either be all possible 
solutions or a set of partial solutions that cover the actual solutions. As stated in 
Theorem 2(2), to solve classical CSP's, fcost should be defined as the constant 
function 6 G 5ft and the parameter o should have the value =. In Table 1 this 
case is given in the first row. 

As well, as shown in Theorem 2(3), a CSP is solved as a COP by instantiating 
o as either > (for maximisation problems) or < (for minimisation problems). The 
value S should be instantiated to the initial cost value from which an optimal 
solution must be found. Traditionally, the range of the cost function (i.e., £)<) 
is the domain 5ft. Rows 2 and 3 in Table 1 show how to initialise both 5 and o in 
these two cases. 



CSP Type 


fcost 




o 


5 


Classical CSP 


constant 






fcost (S) 


Typical Minimisation COP 


any cost function 




< 


T w 


Typical Maximisation COP 


any cost function 




> 


T R 


Max-Min COP 


any cost function 


K x 3? 


< 


(Tir, -Lift) 



Tablel. CSP type depends on parameters instantiation 



In contrast to typical COP's that maintain cither a lower bound or an upper 
bound, our schema also permits a mix of the maximization and minimization 
criteria (or even to give priority to some criteria over others). This is the case 
when D < is a compound domain. Then the ordering in D < determines how the 
COP will be solved. 

Exam-pie 4- Let (V,T>,C) be a COP, D < the domain 5ft 2 with ordering (a, b) < 
(c, d) a < cAb > d, fcost(S) = (fcost 1 (S) , fcost 2 (S)) a cost function on 5ft 2 

for any S € T> where fcost 1 ,fcost 2 T> — > 5ft. Then, if 8 and o are as shown in 
Row 4 of Table f , this COP is solved by minimising fcost ± and maximising fcost 2 . 
However, if < is defined lexicographically, this COP is solved by giving priority to 
the minimisation of fcost 1 over the minimisation of fcost 2 e.g. suppose Si, S 2 and 
S3 are solutions with costs (f. 0,5.0), (3.0, f.0) and (f. 0,8.0), respectively. Then 
with the first ordering the optimal solution is S3 whereas with the lexicographic 
ordering Si is the optimal solution). 

7 Concluding remarks 

This paper analyses the branching process in constraint solving. We have pro- 
vided a generic schema for solving CSP's on finite or continuous domains as well 
on multiple domains. We have proved key properties such as correctness and 
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completeness. We have shown how termination may be guaranteed by means of 
a precision map. We have also shown, by means of an example, how, for systems 
supporting multiple domains, the precision map can be used to normalise the 
heuristic for variable ordering. 

By using a schematic formulation for the branching process, we have indi- 
cated which properties of main procedures involved in branching are responsible 
for the key properties of constraint solving. For optimisation problems, we have 
also shown by means of examples that, in some cases, the methods for solving 
CSP's depend on the ordering of the range of the cost functions. 

By combining a filtering function satisfying our conditions with an appropri- 
ate instance of our schema, we obtain an operational semantics for a constraint 
programming domain (for example: FD, sets of integers, Booleans, multiple do- 
mains, ...,etc) and systems designed for constraint solving such as clp(FD) [?], 
clp(B) [?], DecLic [?], clp(B/FD) [?], CLIP [?], Conjunto [?] or CLP(BNR) [?]. 

Further work is needed to consider how to construct an efficient implemen- 
tation 5 . Moreover, it would be useful to examine how the efficiency of a COP 
solver in our schema could be improved by adding constraints fcost(S)o5 to the 
original set of constraints for solving C, so that exhaustive search is replaced by 
a forward checking mechanism. 

Appendix: Proofs 

A path q G (Natural\{0})* is any finite sequence of (non-zero) natural numbers. 
The empty path is denoted by s, whereas q . i denotes the path obtained by con- 
catenating the sequence formed by the natural number i with the sequence 
of the path q. The length of the sequence q is called the length of the path q. 

Given a tree, we label the nodes by the paths to the nodes. The root node 
is labelled e. If a node with label q has k children, then they are labelled, from 
left to right, q . 1, . . . , q . k. 

Definition 9. (Search tree). Let S be a constraint store for (V,T>,C), e G 9? and 
p G 3?X. The search tree for branch(C, S,p, e) is a tree that has S at the root node 
and, as children, has the search trees for the recursive executions of branch /4 as 
consequence of reaching Line 8 of Figure 1. 

Given a search tree for branch(C, S,p, e), we say that S e = S is the constraint 
store and p e = p the precision at the root node e. Let S q be the constraint 
store and p q the precision at a node q. If q has k > children q . 1, . . . ,q . k, 
then S q is consistent and, if S? = filtering -p(C, S q ), then S? is divisible so that 
choose(S^) = dj and, for some k > 0, split D (dj) = (dji,... ,djk). Then we 
say that S q .i = S^[dj/dji] is the constraint store and p qA — precision(S^) the 
precision at node q . i, for i G {1, ... , fc}. 

5 For example, in COP problems, the top of P can be removed after Line 4 of the 
schema so that P contains only the optimal solution found so far and memory can 
be saved. 
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Lemma 1. Let choose /l be a selecting function for V, split D /l a splitting func- 
tion for D £ Ad, S — (di, . . . , d n ) a consistent and divisible constraint store for 
(V,V,C), dj = choose(S), dj C D and (dji, . . . ,djk) = split D (dj). Then 

(a) Mi £ {1,... ,k} : S^/d^] ^ s S. 

(b) Also, if S' is a solution for (V,T>,C) and S' -< s S, then 

3ie{l,... ,k}:S' ^ s Sidj/dji}. 

Proof. We prove the cases separately. 

Case (a). By Definition 5, #dj > 1 and, by the contractance property shown 
in Definition 6, for all i £ {1, ... , k} dji C dj. Therefore, by Definition 2, for all 
ie{l,... ,k} S[dj/dji] -< s S. 

Case (b). By Definition 5, jfdj > 1 and, by the completeness property of 
the splitting functions shown in Definition 6, 

Vs £ dj . 3i e {1, . . . , k} : s € d jt (1) 

Suppose that S' = (d[, . . . , d' n ). By Definition 2, for alii G {1, . . . , n} d\ — 
{s-} and also, as S' -< s S, s'j £ dj. As consequence, by (1), 3i £ {1, . . . , k} : s'^ £ 
dji so that e?^ C dji. Therefore, by Definition 2, S' <s S[dj/dji\. ■ 

Theorem 1 on page 5. 

Proof. (Property (1). Termination) In the following, we show that the search 
tree for branch{C, S,p,e) is finite so that the procedure effectively terminates. 

Let S e — S and p e — p. If the search tree for branch(C, S e ,p e ,e) has only 
one node then the procedure terminates. Otherwise, the root node e has k chil- 
dren with constraint stores Si where i £ {1,... , k) and S'j = S[[dj / dji]. By 
Lemma 1(a) and Definition 4, for all i £ {1, ... , k}, Si -< s S e and, by Proposi- 
tion 1, precision(Si) <sri precision (S e ). Then, precision(Si) <sri Tsrx- Suppose 
now that precision(Si) = (Tsr, n) for some n £ Integer. Then the test in Line 2 
Pi — precision(Si) < (e, 0) holds and the node containing Si has no children. 
Otherwise, for alH £ {1, . . . , k}, 

Pi - precision(Si) >m (e, 0) (2) 

and there exists some constant I £ 5ft such that 

precision(Si) <srj {(■ x e, 0). 

We show by induction on the length j > 1 of a path q in the search tree that 

precision(Si) — precision (S?) >ri ((j — 1) x e,0). 

It follows that j < I and that, all paths have length < I + 1 (since the second 
condition in Line 3 in Figure 1 holds) and thus there are no infinite branches. 
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The base case when j = 1 follows from (2). Suppose next that j > 1 and 
that the hypothesis holds for a path q of length j — 1. Let q . i q be a child of g 
of length j. Then, by the condition in Line 3 of the if sentence, 

P q .i q - precision(S } q tii ) > m (e,0), 

However, by the inductive hypothesis, 

precision(Si) — precision (5^) >j?x (0 — 2) x e, 0) 

so that, as precision (S*) is 

precision(Si)-precision(S q . iq ) >^x (e, 0) + ((j - 2) x e, 0) = ((j - 1) x e,0).B 

Proof. (Property (2). Completeness) Let i? be a solution for (V, T>, C). Then, 
R ^ s S £ and, by Definition 4, R ^ s 5/. If R = 5/, then 5/ is consistent and 
not divisible by Definition 2, tests in Lines 2-3 hold and R is pushed in the stack 
P. Otherwise, R -< a 5/ and, by Definition 2, 5 1 / is consistent and divisible. As 
Pe = Tsrx, then condition in Line 3 does not hold and node e has k children. 
By Lemma 1(a) and Definition 4, for any q of length m > 1, S 1 ^ f -< s 5^ and 

by Proposition 1, precision(S^) — precision{S q i ) > (0.0) that means that the 
condition p qAq — precision (S q i ) < (e, 0) in Line 3 never holds. It follows that all 
the branches in the tree terminate either with an inconsistent store (because test 
in Line 2 does not hold) or with a consistent and not divisible store (as result of 
holding tests in Lines 2 and 3). Now, we show by induction on the length j > 1 
of a path q in the search tree that 

R <s S{ 3i q e {1, • • • , k} : R ± s Sl lq . (3) 

As, by hypothesis, the procedure terminates then the search tree is finite and 
it follows that there exists some path p = q . q' such that R — S*. Then, by 
Definition 2 5^ is consistent and not divisible and tests in Lines 2 and 3 hold 
so that R is put in the stack P. In the base case, when j = 1, Si = S([dj/dji] 
(i G {1, . . . , k}) and by Lemma 1(b) and Definition 4 3i £ {1, . . . , k} : R ^ s S{ . 
Suppose next that j > 1 and that the hypothesis holds for a path q of length j — 1 
so that R ^ s Sf.IfR ~< s S? then, by Definition 2, S? is divisible so that the node 

S^ has k children by Lemma l(b)and Definition 4 3i q G {1, . . . , k} : R ^ s S qi . 

m 

Proof. (Property (3). Aproximate completeness) Let R be a solution for (V, V, C) 
so that R < s S f _ and, by Definition 4, R < s 5/. Since e > 0.0, and as shown in 
Theorem 1(1), all paths have length < £+1. Therefore, as shown in completeness 
proof, by following (3), there must exists some path q with no children and 
length j > 1 such that R < s S?. If R — S* then R is put in the stack since, by 
Definition 2, R is consistent and not divisible and thus tests in Lines 2 and 3 
hold. Otherwise, as shown in termination proof, the node S? has no more children 
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since the test p q — precision (S?) <sri (e, 0) holds and is put in the stack. 
As R < s S£, by Definition 2, S£ is either a solution for (V, T>, C) or a partial 
solution that covers R. ■ 

Proof. (Property (4). Correctness) Let R be an element in the final state of P. 
As shown in completeness proof, if e = 0.0 the test p q — precision (S^) < (e, 0) 
never holds, for all path q of length m > 1 (also, Line 3 is never satisfied when 
q = e). Therefore, R is in P because there exists a path q where = R is 
consistent and not divisible so that tests in Lines 2 and 3 hold. Moreover, by 
Definition 4, is a solution for (V, T>, C). ■ 

Proof. (Property (5). Approximate correctness or control on the result preci- 
sion) By Theorem 1(4), if R G Po.o then R is a solution for (V,T>,C) and, by 
Theorem 1(2), if R is a solution for (V,T>,C) then R G P .o- Also, by Theo- 
rem 1(3), 3S £l G P ei and 3S £2 G P £2 such that R ds S El and i? <s S £2 . Thus, 
by Definition 3, Po.o dp Pei an d Pa.o dp Pe 2 (also P ei and P £2 are not empty). 

Now we prove that P £l dp P £2 . Suppose that S £2 G P £2 . Then, exists a path 
q of length m > such that S£ = S £2 and S £2 was pushed in the stack because 
tests in Lines 2 and 3 hold so that S £2 is consistent and (a) also not divisible 
(i.e., it is a solution), or (b) the test p q — precision(S^) < (£2,0) holds. 

Suppose (a). If q = e then it is obvious that S £2 is also pushed in P £l . 
Otherwise, for all path q\ with length > 1 and where q = q\ ■ qi the test p qi — 
precision (S^ ) > (e, 0) holds for e — e-i and thus also holds for e = e\. Therefore, 
S £2 is also pushed in P £l . Now suppose (b). Then, if p q — precision(S^) < (ei, 0) 
then S £2 is also pushed in P £l , otherwise since S q covers a solution and as shown 
in proof of Theorem 1(3) there is be some path q' = q . p such that the test 
Pq' — precision{S ql ) < (ei,0) holds and thus S q , is pushed in P £l . By induction 

on the length of p and by (3), it is straightforward to prove that S ql <s S*. 
Thus, by Definition 3, P £l < p P £2 . U 

Theorem 2 on page 7. 

Proof. (Property (1). Termination) This proof is as that of Theorem 1(1). ■ 

Proof. (Property (2)) Observe that if j 'cost (S) = S for all S G T>, then test in 
Line 4* of the extended schema always holds. It is straightforward to prove, in 
this case, that the schemas branch /4 and branch + /4 are equivalent so that all 
properties for the schema branch /4 are also held in the schema branch+ / 'A. ■ 

Proof. (Property (3). Soundness on optimisation) We prove the case when o and 
5 are, respectively, > and -Ld k - The respective case is proved analogously. As 
shown in proof of Theorem 1(2), for e = 0.0, if R is a solution for (V, V, C) then 
there exists some path q of length j > 0, such that R = and the tests in 
Lines 2- 3 hold by Definition 2. Thus, Line 4* is reached for each solution R G 
Sol((V, V, C}), and as consequence, the top of P will contain the first solution 
found that maximises f cost /I. ■ 

Theorem 3 on page 7. 
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Proof. (Property: Approximate soundness) We show that during the execution 
of branch + (C 7 S,p,ei), Line 4* is reached for some S£ top(P £2 ) (where q 
is a path of length mi > 0). As consequence, fcost(S^) d?D < f cost (top (P £2 )). 
It follows that either is in the top of P £l or is not in the top of P £l 
since f cost (top (P £l )) d?D < fcost(S^) so that effectively fcost(top(P £l )) -<d k 
fcost(top(P E2 )). 

Observe that top(P £2 ) is in P £2 because there exist some path q' with length 
m 2 > such that S*, — top(P £2 ) and tests in Lines 2-4* hold. Then, as shown 
in proof of Theorem 1(1), for all path qi such that q' = q 1 . q 2 

p qi - precision (S^) > (£2,0). 

As £i < £2, this also holds in the execution of branch+(C, S, p,Si). Therefore, 
if Sjj, is in P £2 because it is consistent and not divisible, then in the execution 
of branch + (C, S, S", £1) Line 4* is also reached with S^,. Otherwise, S*, is in P £2 

because p q > — precision(S q ,) < (£2,0).. Then, as E\ < £2, and as shown in proof 
of Theorem 1(5), there exists some path q" with length r > mi such that 

p q " — precision (Sy,,) < (£i,0) and S q „ ^ s S ql . 
so that again Line 4* is reached for S*„ . 
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Appendix: A simple example 

Here we show a simple example in the domain Interv (Integer), illustrating the 
flexibility of the schema to solve a CSP in different ways. Let (V, V, C) be a CSP 
where V = {x\,x 2 ,x 3 }, A-p = {Interv (Integer)} and C is the constraint set 

{x\ + x 2 + x 3 < 1, x\ < 1, x 2 < 1, £3 < 1, xi > 0, x 2 > 0, x 2 > 0}. 

Consider also the following cost functions 6 defined on different ranges: 



fcost 3 (xi, x 2 , £3) = (fcost 2 (xi,x 2 ,x 3 ),xi + x 3 ). Range: 5ft 2 
fcost i (xi 1 x 2 ,x 3 ) = (fcost 2 (xi 1 x 2 ,x 3 ),x 2 + x 3 ). Range: 5ft 2 

Consider now the instance branch i n terv(FD)i & s defined in Section 6.1, where 
FD = Integer and assume that choose na ive is as defined in Example 2, p = Tsrx, 
£ = 0.0 and initially the global stack P is empty. Suppose that as filtering algo- 
rithm we define a simple consistency check on the consistency of constraint stores 
in such a way that filtering V (C , S) returns S if S is consistent and the incon- 
sistent store (0, 0, 0) otherwise. Now, assume that branch + (C, S,p, e) is executed 
with different values for 5, o and fcost/1. Since the domain is finite, termination 
is guaranteed even if e = 0.0. Each row in Table 2 corresponds with a different 
execution of the extended schema where 

— Column 1 indicates the way in which the CSP is solved, 

— Column 2 shows the value to which S is initialised, 

— Column 3 the cost function used in the current instance, 

— Column 4 the initialisation of o, 

— Column 5 indicates where is, in the global stack, the solution(s) and 

— Column 6 references the figure that shows the final state of the stack P 7 . 

By simplicity, suppose that during each execution of the extended schema, 
branches are solved by classical backtracking following a classical depth first 
strategy. Then, the CSP is solved in different ways. For instance, to solve the 
problem as a classical CSP (see Row 1 in Table 2), fcost is a constant function 
with value <5 (where S is 1.0) and o is =. Then, all possible solutions for the 
problem are pushed in the stack (see Figure 2(a)). 

Also, Rows 2-3 in Table 2 show how to solve this CSP by maximising and 
minimising the function fcost 2 respectively. The optimal solution is that on the 
top of the stack (see Figures 2(b) and 2(c)). On their turn, Rows 4-7 indicate 

6 The sum is defined to return the mid point of the sum of operand intervals, e.g., the 
sum of two intervals a..b and c.d is (a+c+b+d)/2 that is exactly the mid point of 
the interval a+c..b+d. 

7 To the right of each element in P we write its cost. 



fcost 1 (x 1 ,x 2 ,x 3 ) 
fcost 2 (xi,x 2 ,x 3 ) 



x\ + x 2 + x 3 . 



1.0. 



Range: 5ft 
Range: 5ft 
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CSP Type 


<5 


cost function 


o 


Solution 


Figure 


Classical CSP 


1.0 


fcost x 




Each element in the stack 


2(a) 


Maximisation COP 


J-SR 


fcost 2 


> 


Stack Top 


2(b) 


Minimisation COP 


T K 


fcost 2 


< 


Stack Top 


2(c) 


Max-Min COP (i) 


(J-jf, Tj}) 


fcost 3 


<1 


Stack Top 


2(d) 


Max-Min COP (ii) 


(_Ljf, T») 


fcost 4 


<1 


Stack Top 


2(e) 


Max-Min COP (iii) 


(.Lift, T») 


fcost s 


<2 


Stack Top 


2(f) 


Max-Min COP (iv) 


(J-jf, Tsr) 


fcost 4 


<2 


Stack Top 


2(g) 



Table2. Different solvings of the CSP 



Solution (S) 


fcost^S) 


fcost 2 (S) 


fcost 3 (S) 


fcost 4 (5) 


(1,0,0) 


1.0 


1.0 


(1.0,1.0) 


(1.0,0.0) 


(0,1,0) 


1.0 


1.0 


(1.0,0.0) 


(1.0,1.0) 


(0,0,1) 


1.0 


1.0 


(1.0,1.0) 


(1.0,1.0) 


(0,0,0) 


1.0 


0.0 


(0.0,0.0) 


(0.0,0.0) 



Table3. Evaluation of the solutions to the problems 



how to mix optimisation criteria to solve the CSP. For instance, assume the 
following two ordcrings on 5i 2 : 

( a , b) <i (c, d) a > c and b < d; 

(a, b) <2 (c, d) ^=^> a > c or a — c and b < d. 

Then, row 4 corresponds to the problem of maximising x\ +x 2 +x 3 and minimis- 
ing x\ + xz whereas row 5 corresponds to the problem of maximising x\ + x 2 + x 3 
and minimising x 2 + x 3 . Also, row 6 corresponds to the problem of firstly max- 
imising xi +x 2 +x 3 , and if this cannot be more optimised then minimise x\ +x 3 
(this is consequence of the ordering < 2 ) whereas row 7 does the same but min- 
imising x 2 + x 3 . Figure 2 shows the final state of the global stack for each of 
these cases (also Table 3 shows the evaluation of each solution to the CSP by 
the different cost functions 
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Figure2. The final state of the global stack P in the different solvings of the CSP 



